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Abstract 

A theoretical study of an adsorpti ve reactor which combines tnuitibed pressure swing adsorption 
and chemical reaction is presented; such a reactor is referred to as a pressure swing reactor, or PSR. 
Studies have concentrated on an asymptotic case in which there is the ideal propagation of concen- 
tration waves within the reactor beds; the method of characteristics was employed in the solution of 
the governing PSR equations. The studies assessed the effects of operating conditions, and cycle 
configurations, on the PSR performance. Calculations indicate enhanced reactant conversion when 
compared to conventional steady state plug flow operation. In particular, for some reversible reactions, 
substantial improvements over equilibrium yields have been calculated. For example, for the disso- 
ciation reaction 2A B + C, and where B is the only adsorbing component, approximately two-fold 
improvements o ver the equilibrium yield of product B have been predicted. Such reaction enhancement 
can be attributed to the limitation of the backward traction, which results from the separation of the 
product species B and C 

In addition to the method of characteristics, a cells-in-series method for the asymptotic case has 
been developed, and found to yield calculations consistent with the method of characteristics solutions. 
In a third numerical approach, the spatial discretisation technique of orthogonal collocation on finite 
elements was applied to the governing PSR equations, and the resulting system of ordinary differential 
equations solved by a standard integration algorithm. In this case, many of the simplifying model 
assumptions were relaxed, allowing, for example, the simulation of a non-isothermal PSR with finite 
mass transfer rates. 

One practical significance of reaction enhancement by pressure swing adsorption is a lower tem- 
perature of operation than in a conventional reactor. This would lead to savings in the energy 
requirements of the reactor, and limit the rate and degree of catalyst deactivation due to coke deposition 
or sintering. 
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1. Introduction 



Adsorptive reactors couple the operations of separation and chemical reaction into a 
single process. One of the first examples of such a process is the chromatographic reactor 
(see. for example, refs. [ 1-3] and extensioas of which include real and simulated moving 
bed processes (sec refs. [4-4S] ). Such reactors have advantages over conventional steady 
flow reactors in that enhanced reactant conversion can be achieved. Consider, for example, 
a reactant A flowing through a packed bed, and undergoing a reversible elementary disso- 
ciation reaction 2A«>BH-C. The reversible adsorption of a product component say B, 
results in the separation of B and C in both axial space, and across the solid and gas phases. 
Subsequently* reaction in the backward direction is constrained, shifting the reaction equi- 
librium towards the further production of B and C In the case of industrially important 
endctherraic reactions (e.g. cracking, dehydrogenation, bydrogenolysis), equilibrium lim- 
itations often dictate the necessity of high temperature operation to achieve realistic levels 
of conversion. Reaction enhancement could allow a lower temperature of operation for a 
given conversion, thus giving savings on the process energy requirements, it could also 
make a catalytic process feasible since a lower temperature of operation would reduce the 
degree of catalyst deactivation due to. for example, coke deposition or sintering. However 
major disadvantages of the chromatographic type reactors include non-continuous and 
relatively complex operation, especially in the case of the moving bed processes. Similarly , 
the use of membrane reactors for achieving simultaneous reaction and separation has many 
problems to overcome before its obvious potential can be realised. The scope for research 
for novel adsorptive reactors is therefore apparent. 

This paper presents a theoretical study of an adsorptive reactor which combines multibed 
pressure swing adsorption (PSA ) and chemical reaction. A schematic diagram of such a 
reactor, which is referred to as a pressure swing reactor (or PSR), is given in Fig. 1. 
Typically, the reactor would be packed with a mixture of an active catalyst for the reaction, 
and a selective adsorbent for the adsorption of one or more of the reaction species. Conven- 
tional steps for separation by pressure swing adsorption are employed; see refs. [7,8]. The 
three steps of pressurisation, product release and deprcssurisation constitute the simple 
cycle. In order to achieve a continuous flow of exit gases, two beds are coupled and operated 
in the sequence indicated in f ; ig. I . Other steps can be incorporated into the simple cycle so 
as to enhance product recovery or bed regeneration. For the latter, one such step is the 
countercurrent purge step, and has been considered in this work. A schematic diagram of 
the purge step, and sequencing for a two bed purge cycle, are also shown in Fig. 1 . Important 
advantages of the PSR include a relatively simple operation, based on extensions of an 
existing industrial process and established mechanical technology, and the continuous 
delivery of product A PSR may therefore be suitable for large scale applications. 

fn the following sections, models for a PSR are presented, and calculations illustrating 
the potential of the process given. Emphasis is given to the reversible dissociation type 
reactions, and in particular to 2A*>B+C. Possible practical examples of dissociation 
reactions include the dehydrogenation reactions of ethane to ethylene, and propane to 
propylene. As noted above, conventionai reactors for these are usually operated at high 
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Fig. 1 . Schematic diagram of a PSR showing the basic steps of operation, (b) Sequencing of the steps for two- 
bed processes employing simple and purge cycles. 

temperatures so as to overcome kinetic and, especially, thermodynamic constraints, albeit 
at the expense of high energy requirements. Catalytic processes for these reactions are not 
generally feasible because of the rapid catalyst deactivation due to coking which would 
occur at these temperatures, ll is important to note that unlike the test case reaction scheme 
chosen For this study, the case of a real dehydrogenation is complicated by the effects of 
the pressure dependence of the equilibrium yield. The equilibrium conversion of 2A ~ B + C 
is pressure independent, and thus provides a simplified test case for the analysis of pressure 
swing reaction. However, where there is a net increase in moles due to reaction, a low 
pressure (or lew average pressure) will be desirable for high reaction conversion, even 
though a large pressure ratio is generally desirable for the maximum separation effect. 
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2. Theoretical 



IS 



Difficulties in the simulation of PSR arise from the non-linear nature of the system, and 
the inherent solution discontinuities ot sleep gradients in concentrations. A further com - 
plexity arises from the step-wise (or ramp- wise) changes in pressure at the feed end of the, 
bed, causing periodic reversals in the gas flow direction within the bed. In this work, three 
numerical approaches to the problem have been adopted: the method cf characteristics 
(MOC), and spatial discretisation methods based on cclls-in-series (CIS) and orthogonal 
collocation on finite elements (OCFE). The first two approaches wen; applied to a simple 
model in which, for example, pressure gradients along the length of the bed were neglected, 
reaction species assumed dilute, and linear adsorption isotherms employed. These conditions 
result in a situation in which there is the ideal propagation of a concentration wave fronts, 

1. e., no dispersive effects on the wave fronts. CIS discretisation was carried out using the 
General Process Modelling System (gPROMS) of Barton and Pantelides [9J, currently 
under continued development at Imperial College. OCFE discretisation provided a means 
for the simulation of models which relax many of the usual assumptions employed in 
describing pressure swing processes. Further details of these numerical approaches are now 
given. 

2. /. Method of characteristics 

Principal assumptions for the simple case model are as follows. (1 ) Plug flow of gases. 
(2) Dilute reaclants in an inert and non-adsorbing carrier. (3) An isothermal process. (4) 
No axial pressure gradients. (5) No radial pressure or concentration gradients. (6) Instan- 
taneous local equilibrium between the bulk gas and adsorbed phases. (7) Linear adsorption 
isotherms. (8) Simple power-law reaction kinetics. (9) Beds consisting of a uniform 
mixture of catalyst and adsorbent particles. Given these assumptions, and for the dissociation 
reaction 2A<=> B +C, total and component material balances for a PSR bed can be written 
as 



and 



respectively, where 

1 dP 

/>'=-— 
P dr 



r^UCl-CtCc/Kc) 



(1) 



(2) 



(3) 

(4) 

(5) 
(6) 
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and where s is the local voidage, p is the solid density (kg m~*) t C is the gas phase 
concentration (mol m~-), £, is the forward reaction rate constant (m 6 kg' 1 mol~ 1 s~ 
Af c is the reaction equilibrium constant, m is the adsorption isotherm gradient fmol kg 1 
Pa~* 1 ) , P is pressure (Pa) . r t is the reaction rai e (mol kg 1 s " 1 } > R is the gas constant, ! is 
time (s) , T is temperature (K),u is the superficial gas velocity (m s~ *), y is the gas phase 
mole traction, z is the axial space coordinate ( m) , and the subscript i a component identifier. 
Modified formulations of Eqs. (4) and (5) yield models for other reaction schemes. 

Boundary conditions during the four steps of operation considered in this work are given 
in Table 1. These define the rate of pressurisation or depressurisation (dP/dr), the exit 
velocity of the carrier gas, w L , and the mole fractions of reactant species at z » 0 or z « L. 

The analytical solution of Eq. ( 1 ) can be readily written as 

a = « I/ l»P'(L-2) (7) 

The component balances are of a hyperbolic form, and can therefore be solved by the 
method of characteristics; see, for example Holland and Liapis [ 10], and more specifically 
Morgan and Kirkby [11]. Application of this method to Eq. (2) yields a set of characteristics 
which describe the characteristic directions, i.e., 



(8) 



(9) 



dz u u L +P'(L-z) 
dr'ft** fit 
and the variation of gas composition with lime along these, i.e., 

Eqs. (8) and (9) are suitable for integration by a standard integration algorithm. It is 
important to note however, that since the reaction terra is a function of the compositions of 
all components, then for a given set of characteristics, and say for component j. gas com- 
position values need to be determined for all other components, i+j, at the given z and r. 
For each component this is achieved by interpolation in the z domain, at a given r, 
between adjacent characteristics. In this work, a Runge-Kutta-Merson integration algorithm 
was used, and linear or non-linear interpolation employed at the end of every integration 
time- step. At the cyclic steady state (CSS ), checks on the component material balances for 
the MOC solutions (over a complete cycle) were found to be less than 2% in error. These 
indicate the good numerical accuracy of the model solution. 

Table 1 

Boundary and other conditions for the PSR process 



Processs step 


dP/di 




V/U-0) 




EVessuri&aitao 


Constant >0 


0 


Gomtmxt 




Product release 


0 


Constant >0 






Depress urisation 


Constam<0 


0 


• 




Purge 


0 


Constant <0 


• 


h 



'Determined from the material balances. 

*Based on the average mole fraction of gas collected doling the product release step. 
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2.2. Cells-in-series: gPROMS 



In the CIS discretisation method, the bed is divided into a series of well-mixed cells. The 
well-mixed assumption eliminates spatial gradients within a cell, whilst a spatial gradient 
is maintained from one cell to the next. The governing partial differential equations of the 
PSR model arc discrctised for each cell, and the resulting system of ordinary differential 
equations solved by a standard integration algorithm. Solution accuracy (or convergence) 
generally increases with the number of cells, as does the plug flow condition. 

For cyclic gas flow in a packed bed, a further consideration is required when using CIS, 
i.e., a particular cell can be in one of four general stales at any instant in time. The four cell 
states represent the direction of pressurisation or depressurisation from each end of a cell; 
material balances for these will vary accordingly. For no axial pressure gradients, as in the 
model described in Section 2.1, gas flow throughout the length of the bed, at any instant in 
time, can only be in one direction. Specifically, the pressurisation and product release steps 
are associated with gas flow from the feed to product ends of the bed (say the forward 
direction) , and the countercurrent depressurisation and purge steps with gas flow from the 
product to feed ends of the bed (say the backward direction). For a given cell c of length 
L ct the total material balance for the cell (based on the assumptions listed in Section 2.1) 
can then be written as 

[±(H ! --Hr)//J-P'-0 ( 10 > 

for gas flow in the forward /backward directions. Similarly, the component material balances 
arc given by 



(11) 



(12) 



(13) 



where for gas flow in the forward direction 

s = ( u c - , y,. r - 1 ~ u cytJ tLe 

and for the backward direction 

s = ( - u c . i y lmC + u c y u + { )/L r 

Boundary conditions are as for those listed in Table I . 

Xs mentioned previously, the above ordinary differential equations (Eqs. ( 10) to ( 13)) 
were solved using gPROMS. This provided a concise means for declaring the model 
equations, upon which discrete actions ( at ihc domain boundaries) could be imposed readily 
from one process step to the next Typical CIS simulations involved 100 cells; these, for a 
wide range of operating conditions, were found to be consistent with the MOC solutions. 
i.e./CSS reaction yields within 3% agreement 

2.3. Orthogonal collocation on finite elements 

In the third numerical approach, the spatial discretisation technique of OCFE was applied 
to the partial differential equations describing a PSR. The resulting system of ordinary 
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differential equations were then solved by a standard integration algorithm based on Gear's 
method. For this case, principal model assumptions arc listed below. ( 1 ) Axial! y dispersed 
is plug How of gases. (2) NTon-dilulc rcactants in an inert carrier. (3) A non-isothermal 

n process. (4) Axial pressure distribution as described by Darcy \s Law ox the Ergun equation. 

ie (5) No radial pressure or concentration gradients. (6) Mass transfer limited adsorption as 

al described hy a linear driving force. model. (7) Linear or Langmuir adsorption isotherms. 

0 (8) Simple power-law reaction kinetics. (9) Beds consisting of a uniform mixture of 
catalyst and adsorbent particles. Different forms of the model have also been developed so 

>. as to accommodate both finite rate and instantaneous reactions. The latter involved material 

11 balances on a mass basis, and on an equation describing the reaction equilibrium constraint; 

I; see Vaporciyan and Kadlec [12]. For both forms of the model, the inclusion of axial 

ie .^8" dispersion effects eliminated any solution discontinuities; the steep gradients were then 

in accommodated by low-order Lagrange interpolation polynomials within a finite element 

ys grid. Typically, a discretisation based on five elements, with four internal collocation points 

:d per clement, was found to be computationally precise. Furthermore > when reduced to the 

^e simple case, finite clement calculations were found io be consistent with the MOC and CIS 

ih solutions. 

. ) Further discussions on the application of OCFE to pressure swing problems are given by 
AJpay [13,14]. 

)) . .:.>£; 

es 3. Simulation strategy 

[Tie PSR can be treated as either a separator or reactor, and performance criteria chosen 

1 j for these. For the former, the degree of component separation, or the degree of enrichment 

of a particular component in one stream relative to another, can be assessed. In terras of a 
reactor, predictions of reactant conversions ot product yields arc useful, as arc the produc- 
tivities of die catalyst; see, for example, Lu et al. [15]. Ultimately, comparisons to the 
I ) ..; performance of conventional steady state reactors are required, i.e., to the performance of 

plug flow and stirred tank reactors. For reversible reactions, the conversions attainable by 
• steady state operation are limited by the reaction equilibria. For pressure dependent reaction 

I ) equilibria, further care is needed when comparing PSR to steady state reactor performance. 

For example, for a reversible reaction in a plug flow reactor (PFR) in which there is a net 
increase in moles, a low operating pressure will shift the reaction equilibrium towards 
. product formation, but will also result in a low reactant residence time; hence an optimum 
pressure for maximum conversion may exist. Finally, it is important to note thai even for 
!y pressure independent reaction equilibria, actual reaction rates for the forward and backward 

: a reactions will be functions of pressure, and thus the rates of approach to the reaction 

1S ' equilibria 

In this work, calculations have concentrated on the simple case described above, and the 
method of characteristics employed for model solution. As mentioned previously, the 
reaction scheme 2A <*» B + C was chosen as a test case, i.e., a dissociation reaction in which 
the reaction equilibrium is not a function of pressure. Adsorption of the reactant A and 
cd product C were assumed negligible, leading to the desired separations between B and C, 

try and B and A. Base case design and operating conditions were chosen for the PSR, which 
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Tabic 2 

Ba#e case PSR design and operating conditions for the reaction 2A**B +C. 



Bed length 


1.0 rn 


Gas flow (mol /cycle) ; simple/purge cycles: 


Bed diameter 


0.1 m 




Solid density 


3000 kg/tc 3 


Feed 1.6*9/1.87* 


Total voidagc 


0.75 


Prodnct 0.182/0.182 


Temperature 


473 X 


Purge 0/0.182 


Upper pressure 


10 bar 




Lower pressure 


I bar 


Adsorption isotherm gradients (mol kg 'Pa "*): 


Step times ( s) ; simple/ purge cycles 


A 0 






3 2.5X10-' 


PressorisaXioo 


6/6 


C 0 


Product release 


12/ lb" 




Dcprcssun Ration 


6/6 


Reaction parameters: 


Purge 


0/6 








k 0.005 m 6 kg" 'moT 's" 1 






K c 0.01563 



depict a typical small scale PSA process; see Table 2. The effect of varying an operating 
condition (from the base case) on the CSS performance of the unit was assessed. Tn 
particular, this performance was measured in terms of the totaJ reaction yield of product B 
( Y Q ) , and by the separation attainable between B and C (5 B /c.) , and B and A ( 5 B/ A ) , where 



'uic-ss 



ess 



total amount of B (mol/cycle) collected ail exit stream s 
total amount of A (mol/cycle) supplied as feed 

_ produ ct 



04) 



(15) 



(,T;/ X))e*hau*t uresm 

and where x, and x } are the cycle-average mole fractions of components i and; respectively. 
In Eq. ( 15), the product stream refers to the gas leaving the bed during the product delivery 
step, and the exhaust stream as the combined gases leaving the bed during the ccuntercurrent 
deprcssurisation and purge steps. S|//— 1 indicates no separation of components i and j 
within the PSR, whereas S ifJ > 1 indicates the enrichment of i in the product stream, and 
5,/j < 1 the enrichment of i in the exhaust stream. For the test reaction chesen in this wort 
and for the adsorption of B only, 5 n/A and 5 B/C are expected to be less than unity. Finally, 
the calculated Y a were related to the yield attainable at die reaction equilibrium and, for the 
given design and average feed gas flow rate, by a conventional plug flow reactor. For the 
latter, an operating pressure of 10 bar was chosen (i.e., the upper operating pressure of the 
PSR), this resulting in the attainment of the reaction equilibrium at relatively low reactant 
residence limes. 



4. Results 

In Fig. 2a, PSR, PFR and equilibrium calculations of Y B as a function of the catalyst 
acuvily, * r , arc shown, these for the design and operating conditions listed in Table 2. 
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Regions in which PSR yields exceed those attainable by a PFR are indicated. Furthermore, 
possible improvements over the equilibrium yield are indicated, whereas of course, a PFR 
can only approach the equilibrium yield as k t is increased For the case of no adsorption of 
3 (Fig. 2b), no improvement in the PFR (and subsequently equilibrium) yield was found. 
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Fig. 3. Typical cyclic steady- stale spatial profiles of gas composition during ihe purge $ttp. Arrows indicate 
increasing lirne from the onset of the purge step. 

This suggests thai there is no enhancement of reaction solely due to the periodic operation 
of the reactor. 

In terms cf the mode of operation, a cycle employing a countercurrent purge step is 
shown to be superior to simple cycle operation; sec Fig. 2a and the discussions below. As 
mentioned earlier, the purge step enhances bed regeneration and subsequently the separation 
effect. In Fig. 3 for example, typical CSS spatial profiles of gas composition during the 
purge step arc given, these showing the concentration of product B in the purge gas exit- 
end of the bed (i.e., 2 = 0), and the separation of the product components B and C along 
the length of the bed. This enhanced separation effect due to the purge step is further 
demonstrated by considering the simple and purge cycle configurations employed for Fig. 
2a, but with a non-reacting mixture of A and B as feed. For example, for a feed mixture of 
4% (v/v) A and 0.5% (v/v) B, S BfA for the purge and simple cycles are given as 0.000 
and 0.245, respectively. Hence, greater enrichment of B in the exhaust stream (and, in fact, 
complete separation in this case) when employing a purge cycle. 

Fig. 4a presents calculations of Y B as a function of the important parameters of purge 
amount and the adsorption isotherm gradient of B, ro B . Considerable improvements ovct 
the equilibrium yields are indicated, i.e., approximately 80% improvement over the equi- 
librium yield. For maximum y B , an optimum purge amount is shown to exist For ihc two 
sets of data shown in Fig. 4a, the higher product yields forw B = 2-5x 10~° mol kg" 1 Pa - 1 
correspond to the lower separation factors (thus greater effective separation between com- 
ponents) given in Fig. 4b. The optimum purge amount, however, does not correspond to 
any minimum in the separation factor-purge amount plots. For the case of no adsorption of 
B (calculations not shown), Y H falls below that attainable at the reaction equilibrium and 
by PFR operation; S B/A and Sb, c equal unity for this case. 
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Fig. 4. (a) Total yield of product B and (b) separation factors as a function of the purse 
ter°ffloX10 7 . 



amount; para me- 



ld Fig. 5. calculations of Y n as a function of the combined effects of the purge amount 
and the cycle time settings are shown. In this case, the different cycle time settings corre- 
spond to different feed gas and purge gas flow rates. For example, for a given purge amount, 
any increase in the duration in the purge step will result in a decrease in the purge gas flow 
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Fig. 5. Tola! yield of product B as a function of the purge amount; parameter^ process step time seeing (s) for 
the prcssurisation, product delivery, deprcssurisation and purge steps, corresportdb^y. 

rate during this step. The results indicate an optimum purge amount, this a function of the 
product gas flow rate. 

The relationship of the separation factor to reaction conversion is farther demonstrated 
by considering the effect of the product amount on Y B , this again for a purge cycle. In Figs. 
6a and 6h Y B is shown to decrease with increasing product amount, whilst S B/C and 
5 B/A approach unity. This loss in separation performance with increasing product amount 
is in agreement with conventional PSA operation; see, for example, Kirkby I i6] and 
Espatalier-Noel [17]. 



5. Discussion 

The enhancement of reaction by PSA has been indicated by the above calculations. For 
a reversible dissociation reaction, separation of the product species along the length of a 
PSR bed leads to a reduction in the backward reaction rate ( e.g. B + C -* 2A ) , and subse- 
quently enhancement of the net reaction in the forward direction. The inclusion of a purge 
step is shown to be particularly beneficial in enhancing this separation, and hence reaction. 
Whilst component separation factors provide a useful indication of the degree of reaction 
enhancement, the complex relationships between the gas residence time and the degree of 
separation and reaction precludes any quantitative relationship. This was demonstrated in 
Figs. 4a and 4b in which optimum conditions could not be correlated to the separation factor 
data. 

In addition to the spatial separation of the reaction species, a second possible mechanism 
responsible for reaction enhancement can be attributed to the transfer of a product species, 
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Rg. 6. (a) Total yield of product B and (b) separation factors as a function of die product amount for a purge 
cycle. 

say B, from the gas phase to the adsorbed phase through adsorption. The depletion of B 
from the gas phase results in further reaction for product formation. During the depressur- 
isation step, the relatively large gas flow rates, and hence short gas residence times for 
reaction, then allows the collection of desorbed B in the exhaust stream. Preliminary 
calculations in this work have demonstrated this effect for the single product reaction A ~ B, 
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in which the adsorption of B was again found to give PSR yields greater than that attainable 
at the reaction equilibrium, i.e., yields of approximately 5% in excess of the equilibrium 
yield. 

As mentioned above, real dchydrogenation reactions are complicated by the effects of 
pressure on the reaction equilibrium. Operating regions in which a PSR exceeds the best 
steady state PFR yield ( for a given design) . or highest equilibrium yield, are currently being 
sought for such reactions. As experimental verification for the enhancement of catalytic 
reaction by PSA. a multimodal lab-scale PSR is also being constructed. The dehydrogenation 
of methylcyclohexanc to toluene has been chosen as a test reaction in the first instance. A 
commercial Pt -AI 2 0 3 catalyst will be used and admixed with a selective zeolite or siHcalite 
adsorbent. For this real reaction and adsorption system, any effects due to, for example, 
pressure gradients or mass transfer limitations, will be assessed using the OCFE based 
model described in Section 2.3. 

Finally, theoretical studies will continue to assess die effects of operating conditions and 
modes of operation on the PSR performance. Explorative calculations for the optimum 
operating conditions for a purge cycle PSR are also being carried out. Further potential 
areas of research include pressure swing reaction for enhancing product selectivity, and 
with the advent of zeolitic catalysts, the use of these combined adsorbents/ catalysts in a 
PSR. 



6. Conclusions 

Models of a novel reactor which combines separation by PSA and catalytic reaction have 
been developed, and tested for their numerical precision. For an asymptotic case, in which 
there is the ideal propagation of concentration wave fronts, and for the reversible dissociation 
reaction 2A <=> B + C. almost two-fold improvements over the equilibrium yield of B have 
been predicted. When compared to conventional steady flow reactors therefore, potential 
advantages of a PSR include a higher reactor productivity, or a lower temperature of 
operation for a given conversion. The latter will lead to savings in the energy requirements 
of the process, and prolong the activity of the catalyst or, indeed, could make a catalytic 
process feasible in cases where conventional catalytic reactor operation would be impossi- 
ble. When compared to early examples of adsorptivc reactors, e.g. the chromatographic 
reactors, a PSR offers process simplicity and the continuous delivery of product gas, and is 
thus suitable for large scale operation. Experimental verification of reaction enhancement 
by a PSR is currently underway: the more facile dehydrogenadon reaction of melhylcy- 
clohexane to toluene ( atca. 473K) has been chosen in the first instance, for which a mixture 
of a selective zeolite adsorbent for toluene and a platinum based reforming catalyst is being 
employed. 



Notation 

C gas phase concentration, mol m " 3 

kf forward reaction rate constant m* kg"" 1 mol * * 1 s " 1 
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L 
K 

m 
P 
P' 

'x 

R 

5 

t 

T 

u 

A 



reaction equilibrium constant 
bed length, m 
cell length, m 

adsorption iscthenn gradient, mo! kg" 1 Pa" 1 
pressure. Pa 

pressure parameter, s ~~ 1 (seeHq. (3)) 
reaction rate, mol kg " 1 s~ 1 
gas constant, = 8.314 J mol 1 K 1 

cell state equation for concentration gradient; sec Eos. (12) and (13) 
separation factor for components i andj; sec Eqs. (15) 
time, s 

temperature. K 

superficial gas velocity, m s" 1 
average mole fraction 
y gas phase mole fraction 
Y total reaction yield; see Eq. ( 14 ) 
z axial space coordinate, m 
Greek 

P adsorption parameter ( see Eq. (4 ) ) 
e total bed voidage 
v power law coefficient 
p solid density, kg m " 3 
Subscripts 

0, L feed end condition, product end condition 

A, B, C, i component identifiers 
c ceJI identifier 
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